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Abstract 

A statistical thermodynamic approach of moving particles forming an elastic body is presented 
which leads to reveal molecular-mechanical properties of classical and nonextensive dynamical 
systems. We derive the Boltzmann-Gibbs (BG) entropy form and relate it to Newton's law of 
motion in relation to a distinct tensile force acting on the systems at constant volume and number 
of particles. Tsallis generalization of the BG entropy is deduced assuming the thermal energy of 
the particles to be proportional to their energy states by the nonextensivity factor q — 1. 
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I. OVERVIEW 



It was pointed out long ago by Einstein that the foundations of statistical mechanics 
should lie on the underlying dynamics All statistical mechanics results should derive 
from the dynamics if the systems considered are systems consisting of moving particles. 

Within the classical Boltzmann-Gibbs (BG) statistics, it is still unkown how to derive 
rom say Newton's law of motion the remarkable expression for the thermodynamic entropy 
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S B G/k B = - J^Pi In Pi = ( hi — ) , (1) 

U. \ Pi/ 

where mean values < (• • •) > correspond to J2i{' ' ') Pi- Even if this functional form has 

been profusely verified on isolated systems in thermal equilibrium, a rigorous derivation on 

how Sbg descends from dynamics has been limited to a dilute gas only jlj]. It is usually 

approximated via the statistics of a large ensemble of identical replicas (see, e.g., [^j]). Only 

recently, there has been computational evidence on the fact that BG thermal equilibrium 

distribution descends from Newton law j^j]. 

Also unclear at present , is the connection between the dynamics of a particle system 

and the Tsallis nonextensive generalization of the BG entropy 

/ 1 \ w 1 1 - T w r> q 

SJk B = (in, - )=5> In, - = ^= lft , (2) 

\ Pi/ f^i Pi 9-1 

where 

1-q 1 

, _ Pi ~ 1 /oN 

In? Pi = 1 _ , (3) 

which for q = 1, it coincides with Sbg giving positive Lyapunov exponents p. The nonexten- 
sivity is given by the pseudo-additivity property S q (A+B) = S g (A)+S g (B) + j^ L S q (A)S g (B); 
A and B being two independent systems. 

In order to shed some light on this physics dilemma, we propose a simple statistical 
thermodynamic model where moving particles of known type are assumed to form an elastic 
specimen (or body). It is shown that this type of approach leads to reveal intrinsic molecular- 
mechanical properties on classical and nonextensive dynamical systems in relation to a 
distinct tensile force acting on these systems at constant volume and number of particles. 
In particular we derive the entropy of Eq. (pQ) and relate it to Newton's second law F = ma. 
A new scenario for the entropic index q in Tsallis statistics in terms of the energy of the 
system is proposed. 



II. FORMALISM 



Let us consider a system of non-interacting moving particles and let pi be the probability 
that the system is in the microstate i(= 1, • ■ ■ , W) given by the product of two positive 
functions Ui and Wi and satisfying the normalization condition 

w w 

^Pi = ^2uiWi = 1 . (4) 
i=i i=i 



Derivative of this relation with respect to changes in particle displacements x, considering 

|| and acceleration a = ||. 



particles of mass m to move at velocity v = % and acceleration a = gives 



f ± dx ~[ dx 

or, alternatively, 
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f^ur'F^w, ■ (6) 

Summing up both sides by J2i=i(^r) In u% it is straightforward to show that 
d dv ■ 

_<ln^>=--<lnt*>+£(^)lnp, . (7) 

Let us also consider the purely thermodynamic relation for the Helmholtz free energy 
A = E — TS, assuming the temperature T and volume V to be independent variables 
with E the internal energy, A the Helmholtz free energy and S the entropy of the system. 
Therefore from Eq.fJZj) we can readily identify (or, alternatively, define) three quantities 

w 

A/k B T -> < ln Ui >= ^Pi ln u { , (8) 

i=l 

W 

E/k B T — > - < ln Wi >= ~^Pi In Wt , (9) 

i=l 

with A;^ a constant added to have reduced units. The choice of the sign for E is to ap- 
proximate the system energy as the sum of individual contributions (or mean energy) 
E ^< E >= J2iPi e i, where Ci/ksT = — ln Wi represents their energy states, and thus 
obtain 

Wi = exp(-ei/k B T) >0 . (11) 
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This relation is also found to be the solution of the differential equation: " "V = — 1. 
Furthermore from Eqs. (JSJ) and (fTU|) it follows immediately that the entropy satisfies 

= • (12) 

which is the basic result stated by the BG entropy of Eq.flTJ). 

Since the moving particles are assumed to form an elastic body, then the total work along 
an infinitesimal elongation dt is given by dW = dA = fp di for a change at constant T, with 
/ a (compressive < 0, or expansive > 0) tensile force in the direction of I. This relation 
is obtainable directly from the first and second laws of thermodynamics in a reversible 
isothermal process p, assuming that there is not another force acting on the system in 
addition to ft as for example due to hydrostatic pressure P (dW = —PdV ~ 0; recall 
volume V here is an independent variable). Therefore, the dynamics of our particle system 
relating thermodynamic quantities reduces to 

"(§a-^i(g»H-(f),- r (i),- <i3) 

Thus through Eqs.(jHJ) and (jHJ) the sytem entropy can be correlated to Newton's law of motion 
in relation to a distinct tensile force acting on the system. 



III. REMARKS 

At a first glance the choice for pi in Eq. (pij may be seen quite arbitrary, but in view of the 
results derived for S, A and E starting only from this quantity it seems to be reasonable. As 
shown below, this choice can further be justified in relation with the measure for an inverse 
temperature for the system. Using Eqs.(JH), 0, (fTTj) and (fT2j). we obtain at once 

~ 1 + TTT— = 1 - — ln u w , (14) 



A(E/k B T) ln w w x . „ , 

where A represents functions difference between consecutive microstates fluctuations. 

Although we do not know how each of our and Wi functions may depend explicitly on 
the displacement x, we can still estimate all thermodynamic relations from given average 
values of either u or w at each state i. For example, interesting to note are the results 
obtained when setting m = 1/Z, a sort of mean value for u (Vz) and Wi variable. From Eqs.(jl} 
and (fTTj) . it follows then that pi = 6 £ ^ gT , which corresponds to the Gibbs distribution in 



thermal equilibrium and satisfies ^ = Using the normalization condition for pi 

we derive in turn the partition function Z = J2i e~ € ^ kBT , where the summation is over all 
microstates of the system. From this result and that of the energy states q plus Eq. ([§]). we 
obtain a relation for the mean energy of a system at temperature T as < E >= — ( g^n^m ) • 
For the Helmholtz free energy of Eq.fjHJ) we can write A/k B T = — In Z, which is also a well- 
known classical statistical mechanics result 

Let us analyse next a second case, i.e., the case of having an averaged w (Vz) and it, as 
the variable. To derive Tsallis statistics we assume the thermal energy of the particles k B T 
to be proportional to their energy states e by the nonextensivity (integer) factor q — 1 for 
all the z'-microstates. From Eq. ljllj) and q ^ 1 this is equivalent to consider 

In Wi = . (15) 



Using the ln 9 definition of Eq.(j2J) in conjunction with Eqs.(jU) and (fTo} . we then get the 



relation — ln g (1/pi) — In (pi/ui) = ~^z\ - For q — > 1 and by Taylor expansion [2] we can also 
approximate p\ = PipJ^ 1 ^ = Piexpihi p^ 1 ^ 9 ^) ~ Pi[l + (1 — q) In pj 1 ]. Thus in this limit 
ln g p~ x « In pj . These results imply that 

In u % « ^— j- . (16) 
Therefore within our statistical thermodynamic approach it follows that 

A/k B T^J2Pi{^Zi) ' ^ T ~-f>(y^) = e A* r , (17) 

whence 

S/k B = E/k B T - A/k B T n ^ Pi , (18) 

which corresponds to the entropy term of Eq.(|2J) introduced by Tsallis. As in Ref.Q], if we 
assume equiprobability Pi = 1/W (Vz) we also obtain S = ln g W. 

Possible values of the q index are constrained by the logarithm functions of Eqs. (jl5|) and 
(J16|) . The condition < Pi = UiWi < 1 (Vz) implies that < exp(l/(l — q)) exp(pf^ 1 / (q — 
1)) < 1. Hence we can deduce that @i(q) = {p q l — l)/(q — 1) < 0, due to the positive and 
continue nature of the exponential function. Then the limiting behaviour of this function 
as q — > 1 is obtained by L'Hopital's rule because this is an indeterminate form of type 
— . We found lirn g _>i ©j(g) = hipi, hence < exp( Inpi ) < 1. From this result we also 
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have < pf^ 1 < 1 (Vi, Vg), which via Eq. (J15|) it implies that hiWi < &i(q) < 0. This, in 
turn, implies a positive Wi in agreement with Eq. (fTT|) such that < Wi < 1. Therefore we 
can conclude that the range of validity of our statistical thermodynamic model for Tsallis 
statistics is q e (±oo, 1]. 

On the other hand, using Eqs. lfTo]) and (|H)j) for (u,w) at the upper microstate i = W, 
Eq. (jl4j) for q ^ 1 takes the form 



Pi" 1 • (19) 



A(£/fc s T) 

This implies that the thermodynamic definition for the temperature at constant volume and 
number of particles T = follows for g >> 1 since, in physical terms, p w < 1 for 

large subextensive systems. For superextensive systems with integers q < 1, this relation 

A(5/fc B ) 
A(£/fc s T) 



becomes negative. For q — 1, from Eq.((TU) we readily found xt^ttt] = 1 as should be 



expected since 5*^=1 recovers Sbg- 

Finally for the force on the system within Tsallis statistics we have 

^■^(S 1 )^™ • (2o) 

which gives a physical significance to the g-values representing complex systems and corre- 
lates them to the system dynamics. 

In summary, the new insights gained by our statistical thermodynamic approach of mov- 
ing particles forming an elastic body are on the molecular-mechanical properties derived 
for classical and nonextensive dynamical systems including a possible correlation between 
force and entropy in the system. Our initial factorization for pi = UiWi may still be seen 
as rather arbitrary, but we have shown that this simple multiplicative form leads to in- 
teresting connections between an applied stretching force (tension) and thermodynamics 
quantities of dynamical systems []]]. Such class of normalized product of positive functions 
appears formally, e.g. in the analysis of stochastic processes on graphs according to the 
Hammersley-Clifford Theorem. 
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